Astronomy Sz Astrophysics manuscript no. aa2872 


February 2, 2008 


(DOI: will be inserted by hand later) 





in 
o 
o 

(N 

< 



> 

o 

o 
in 
o 

S3 

Oh 
6 



Very compact radio emission from high-mass protostars. 
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Abstract. This paper reports 43 GHz imaging of the high-mass protostars W 33A, AFGL 2591 and NGC 7538 IRS9 
at ~ 0.04" and ~ 0.6" resolution. In each case, weak (~mJy- level), compact (0~1OO AU) emission is detected, 
which has an elongated shape (axis ratio ~3). However, for AFGL 2591 and NGC 7538 IRS9, the emission is 
single-peaked, while for the highest-luminosity source, W 33A, a 'mini-cluster' of three sources is detected. The 
derived sizes, flux densities, and broad-band radio spectra of the sources support recent models where the initial 
expansion of H II regions around very young O-type stars is prevented by stellar gravity. In these models, accretion 
flows onto high-mass stars originate in large-scale molecular envelopes and become ionized close to the star. These 
models reproduce our observations of ionized gas as well as the structure of the molecular envelopes of these 
sources on 10 3 -10 4 AU scales derived previously from single-dish sub-millimeter continuum and line mapping. 
For AFGL 2591, the 43 GHz flux density is also consistent with dust emission from a disk recently seen in near- 
infrared 'speckle' images. However, the alignment of the 43 GHz emission with the large-scale molecular outflow 
argues against an origin in a disk for AFGL 2591 and NGC 7538 IRS9. In contrast, the outflow from W 33A 
does not appear to be collimated. Together with previously presented case studies of W 3 IRS5 and AFGL 2136, 
our results indicate that the formation of stars and stellar clusters with luminosities up to ~10 5 Lq proceeds 
through accretion and produces collimated outflows as in the solar-type case, with the 'additional feature' that 
the accretion flow becomes ionized close to the star. Above ~10 5 Lq, clusters of H II regions appear, and outflows 
are no longer collimated, possibly as the result of mergers of protostars or pre-stellar cores. 

Key words. Stars: Circumstellar matter; Stars: formation; Instrumentation: high angular resolution 



• i-h , 1. Introduction pockets of ionized gas have developed, but stay confined 

/\ • . to the stellar vicinity. (5) Compact, followed by classic H II 

Lj ■ High-mass stars 10 M fi play a major role in the evo- • i ,i • • j j u j j • n 

*3 , ° y J regions, where tnc ionized gas expands nydrodynamically 

. CG , lution of their host galaxies, through strong ultraviolet and disrupts the pa re ntal molecular cloud (for a review, 
radiation and winds, and through supernova explosions. see | Garav fc Li 
However, the origin of high-mass stars is not well un- 
derstood, due to the large distances ( Z kpc) and heavy Despite this progress, key ques tions remain unan- 
extinctions (Ay > 100) of the regions they form in. The swer ed, which hav e been reviewed bv lChurchweil fc002al) 
optically visible main sequence life of OB-type stars is and IVan der Takl (|2003|) . The paramount question is 
preceded by a n 'embedded' phas e which lasts «15% of whether high-mass stars form via (disk) accretion like 
their lifetime l|Churchwel]| l2002~b(l . Observations at mid- low-mass stars, or that other mechanisms such as co- 
infrared through radio wavelengths have shown that this agulation of lower-mass stars or protostellar cores must 
embedded phase can be subdivided into several groups be invoked l|Bonnell et alJll9 98i. The ubiquity of mas- 
of objects: (1) Infrared dark clouds, where internal den- sive bipolar outflows, detected in 2 1 out of a sample o f 
sity maxima and temperature minima likely represent the 26 high-mass protostellar objects bv lBeuther et all (2002) 
initial conditions of high-mass star formation. (2) High- provides strong circumstantial evidence for disk accretion, 
mass protostellar objects, where the central star is sur- The second central question of high-mass star formation 
rounded by a massive envelope with a centrally peaked is the relation with clustered star formation, i.e., the rela- 
temperature distribution. (3) Hot molecular cores, with tion between the mass distri bution of a star- forming region 
large masses of warm and dense molecular gas, and large and its stellar density (e.g.. lTesti et alJll999|) . Addressing 
abundances of complex organic molecules evaporated off these questions requires sub-arcsecond resolution, which 
dust grains. (4) Ultracompact H II regions, where large at sub-millimeter wavelengths, with the Sub-Millimeter 
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Array (SMA), is just coming within reach l|Beuther et alJ 
However, sub-arcsecond resolution can already be 
achieved at centimeter wavelengths, where two important 
diagnostic tools for high-mass star formation are contin- 
uum emission from dense ionized gas and H 2 O maser emis- 
sion. 

Although the relation between the weak radio contin- 
uum emission detected in some high-mass protostars to 
their overall luminosity is presently not understood, it is 
clear that it marks a signpost for the exact position of the 
object. For example, in the case of Orion-KL, the weak ra- 
dio source "I" almost certainly marks the position of the 
exciting object in the region, as the excitation of the sur- 
rounding SiO masers requires e xtreme temperatures and 
densities ijMenten fc Reidlll995|) . Apart from the signpost 
effect, the detection of radio emission, which in the case 
of Orion-I is optically thick at frequencies up to at least 
345 GHz, can also put constraints on the nature o f the 
embe dded protostar, in particular the accretion rate ijTanl 
120081) . 

Maser emission in the H2O 6i6-5 2 3 line at 22.235 GHz 
traces molecular gas with temperatures of ~500 K and 
n(H 2 )-10 9 cm" 3 JNeufeld fc Melnicklll99ll) as occur for 
instance in shocks. The high brightness of the maser emis- 
sion permits high angular resolution techniques and makes 
it a powerful tool for studying gas motions on small scales. 

We are exploring the diagnostic power of high- 
resolution cm-wave observations to study the forma- 
tion of high-mass stars. In the source AFGL 2136 
ijMenten fc van der Takll2004f) . weak, compact radio con- 
tinuum emission was detected, as well as H 2 maser 
emission which is redshifted from the systemic veloc- 
ity by «4 km s . This emission may arise in clouds 
of infalling gas that are heated in accretion shocks. 
In another dee ply embedded high-m ass star-forming; re- 
gion, W 3 IRS5. IClaussen et all l|l994|) and lTieftrunk et all 
( 1997) detected a cluster of seven compact radio con- 
tinuum sources, some of which are transient. Combining 
high-resolution millim eter- wave and mid-infrare d observa- 
tions of these sources. TVan der Tak et all l)2005|) conclude 
that some of them are probably evaporating and/or shock- 
ionized clumps in the ambient material, while others, hav- 
ing mid- infrared counterparts, most likely contain young 
high-mass stars. 

This paper presents new radio observations at high 
angular resolution of three deeply embedded regions of 
high-mass star formation. The sources are a subset of a 
sample of bright mid-infrared sources studied from the 
grou nd llWillner et alJl98llMitchell et alll99(]|) and with 
ISO ( Van Dishoecklll998l) . Single-dish sub-millimeter con- 
tinuum and molecular line mapping bv lVan der Tak et alJ 
( 200fJ showed massive envelopes with a centrally concen- 
trated temperature and density structure. For the present 
radio continuum study, we selected four sources which 
cover a range of evolutionary states within the embedded 
phase, a s indicated by their sub -millimeter and infrared 
spectra l|Va.u der Tak et al .1120031) . 



W 33 A is a highly luminous object (L=lxl0 5 L Q ) 
in the Galactic molecular ring. The kinematic distance 
based on CS and C 34 S lines is 4 kpc l|Van der Tak et 
2000). Its massive, cold envelope («1100 M Q within 
74,000 AU radius) is a favourite target for observations 
of m id-infrared absorp tion features due to icy grain man- 
tles llGibb et alJl200fi). The millimeter emission was re- 
solved bv IVan der Tak et all |2000) into a double source 
with separation 5". SMA observations will be presented by 
Shirley et al (in prep). Two H 2 maser spots are known, 
whose centroid lies 0'.'75 fr om the millimeter conti nuum 
peak at position ang le -49° l|Forster fc Caswellll999j) . The 
maser velocities are redshifted by 0.2 - 1.5 km s _1 relative 
to the dense core. 

AFGL 2591 is a relatively isolated luminous 
(L=2xl0 4 L Q ) protostar in the Cygnus X region. 
Distances in this reg ion range from 0.5 to 2 kp c; we adopt 
d=l kpc and refer to I Van der Tak et alJ l)l999|) for further 
discussion. The source is kno wn for its very high v elocity 
CO mid-infrared absorption (Mitc hell et "all 11 9891) . With 
an envelope mass of only ss40 M Q within a 30,000 AU 
radius, it is probably the least embedded (most evolved) 
source studied here. A combination of compact radio con- 
tinuum emission and H9O masers similar t o AFGL 2136 
has been observed bv iTrinidad et af] (j2Q03). However, in 
AFGL 2591, the maser spots trace a helical structure in 
position-velocity space, probably due to a precessing out- 
flow. 

NGC 7538 IRS9 is a luminous (L=4xl0 4 L Q ) source 
in the NGC 7538 region, which optical spectrop hotom- 
etry puts at <i=2.8 kpc ijCramoton et al 1 1l978h . This 
source combines strong mid-infrared absorption features 
by icy grain mantles with high velocity CO emission 
l)Schutte et alll996tlMi"tchell fc Hasegawalll99l|) . The en- 
velope mass of «430 M Q inside a 66,000 AU radius makes 
it an intermediate object betwe en W 33A and AFGL 2591 
in terms of envelope evolution. iKameva et al ] dl990l) de- 
tected two H 2 maser spots in NGC 7538 IRS9, which 
have their centroid 0'.'54 separated from the millimeter 
continuum peak at position angle 40° . Their velocities are 
blueshifted by 4.7-17.1 km s" 1 from the cloud velocity 
and they probably arise in the outflow. 



2. Observations 

Initial radio observations of W 33A, AFGL 2136, AFGL 
2591, and NGC 7538 IRS9 were obtained on 2001 
September 9, with the NRAO Very Large Array (VLA) 1 
in C-configuration. Twenty-six antennas observed the 
sources sequentially in standard continuum mode: two in- 
termediate frequency (IF) pairs of width 50 MHz each 
(43 MHz 'effective' bandwidth) separated by 50 MHz. The 
central frequency of 43.3 GHz is called Q-band by radio 
engineers. Calibration of the antenna gains and phases 



1 The National Radio Astronomy Observatory (NRAO) is 
operated by Associated Universities, Inc., under a cooperative 
agreement with the U.S. National Science Foundation 
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Fig. 1. VLA C-array 43 GHz images of our targets. Fig. 2. VLA A-array 43 GHz images of W 33A (top), 

Contours start at 0.6 mJy/beam and increase by 0.6 AFGL 2591 (middle) and NGC 7538 IRS9 (bottom). 

mJy/beam for W 33A (top); for AFGL 2591 (middle), Contours start at 0.45 mJy/beam and increase by 0.15 

contours start at 0.7 and increase by 0.7 mJy/beam; for mJy/beam. The boxes in the W 33A image are enlarged 

NGC 7538 IRS9 (bottom), contours start at 0.5 and in- in the next figure, 
crease by 1.0 mJy/beam. 



was done by 'fast switching' to the sources 1817-254, 
1829-106, 2013+370, and 2248+439 every 140 seconds. 
The antenna pointing was checked every hour with a 1— 
minute integration on the phase calibrators at a frequency 
of 8.4 GHz Ox-band'). 

Follow-up observations were made on 2002 March 23 
and 25, with the VLA in A-configuration. Twenty-three 
and twenty-four antennas produced useful data on these 



respective days. The observing procedure was the same as 
before, except that 1829-106 was used as phase calibrator 
for both W 33 A and AFGL 2136. 

The data were edited, calibrated and imaged with 
NRAO's Astronomical Image Processing System (AIPS). 
At the time of these observations, corrections for elevation- 
dependent antenna gain and atmospheric opacity were not 
applied automatically to VLA data, so we applied them 
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Fig. 3. Close-ups on the three 43 GHz sources in our VLA 
A-array image of W 33A. Contour levels are as in previous 
figure. 



ourselves. In addition, small corrections to the positions 
of some antennas were made. 

The absolute flux density scale follows from obser- 
vations of 1328+307 (3C286) which has a flu x density 
of 1.46 Jy at this frequency (jOtt et alJ fl99^ . For the 
phase calibrators we obtain flux densities of 0.49±0.01 Jy 
for 1817-254, 0.63±0.01 Jy for 1829-106, 2.70±0.06 Jy 
for 2013+370 and 0.28+0.01 Jy for 2248+439. In A- 
configuration, the primary flux calibrator, 3C286, is re- 
solved, so we did not use a point source model, but a de- 
tailed model (in the form of Clean components) supplied 
by the VLA staff. We find flux densities of 0.71+0.02 Jy 



for 1829-106, 3.61+0.12 Jy for 2013+370 and 0.29+0.01 Jy 
for 2248+439. 

Figure ^ shows the C-array maps, obtained by a 
Fourier transform of the uv data with uniform weighting 
and deconvolution with the Clean algorithm. These maps 
have rms noise levels of 0.23 mJy/beam. The synthesized 
beam sizes (FWHM) are 0.84x0.41 arcsec at position an- 
gle (PA) -20° for W 33A, 0.53x0.41 arcsec at PA -44° for 
AFGL 2591, and 0.55x0.42 arcsec at PA -19° for NGC 
7538 IRS9. Figure |2 shows the A-array maps, which have 
rms noise levels of 0.18 mJy/beam. Restoring beam axes 
are 78x36 mas at PA -28° for W 33A, 43x37 mas at 
PA -36° for AFGL 2591, and 47x36 mas at PA -42° for 
NGC 7538 IRS9. Note that the regions shown in Fig. El 
are about the size of the compact emission sources from 
Fig. ID 

3. Results 

3.1. Positions and sizes 

Most images show single point sources. Only the C-array 
map of AFGL 2591 has an additional nearby extended 
H II region. This source was seen before by ICampbelll 
l|l984|) at 4.9 GHz and bv IVan der Tak et alJ l|l999|) at 
86 GHz, and appears partially resolved in our data. The 
flux density in a box around this source is 67.6 mJy; 
the source is presumably resolved out in A-array. We do 
not detect the somew hat extended source VLA2, which 
iTrinidad et alJ |2003) only detect at frequencies below 
10 GHz. This source may be of non-thermal nature, or 
resolved out. 

Table ^ presents the results of Gaussian fits to the 
observed compact radio emission, obtained with the task 
JMFIT in AIPS. The A-array map of W 33A may show 
two weak sources in addition to the source seen with the 
C-array. We call the 43 GHz sources in W 33 A Ql, Q2 
and Q3, in order of decreasing brightness. Figure shows 
close-ups of these sources. If sources Q2 and Q3 are con- 
firmed by more sensitive observations, W 33A may resem- 
ble the W3 IRS5 region, where a cluster of radio sources 
exist s, some of which are inter nally, others externally ion- 
ized l|Van der Tak et alJl2005(l . 

The position of W 33 A Ql coincides wi th the 86 GHz 
source MM1 from I Van der Tak et all l|2000l) . The A-array 
image also shows an extension of the emission toward 
the southwest, which was also seen at 86 GHz and called 
MM2. Our data do not clearly s how a second peak, how- 
ever. The H2O maser centroid jForster fc Caswell 1 1999j) 
lies 0'.'44 from the 43 GHz peak, at position angle -78° 
(i.e., due West). The positions of AFGL 2591 and NGC 
7538 IRS9 are consistent with the millimeter-wavelength 
(van der Tak et al. 1999, 2000) and mid-infrared positions. 
The positions found in A-array are consistent with those 
found in C-array within their 2a errors. 

The fitted FWHM sizes of the sources are in columns 
6 and 7, and the orientation in column 8 of Table Q 
The given sizes are deconvolved values, even though most 
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Table 1. Gaussian fits to the observed radio emission. Numbers in brackets are uncertainties in units of the last 
decimal. 



Source 




[divou) 


r^eaK i v 


lotai Dv 


Major axis 


Minor axis 


PA 

m 




hh mm ss 


o / // 


mjy/beam 


mjy 


mas 


mas 


deg 






C- 


-array 










W 33A 


18 11 44.2121(12) 


-17 52 57.793(29) 


3.0(2) 


4.3(5) 


530(140) 


270(110) 


168(31) 


AFGL 2591 


20 27 35.9523(22) 


+40 01 14.694(17) 


2.6(2) 


3.3(5) 


440(110) 


<160 


89(15) 


NGC 7538 IRS9 


23 11 52.8696(14) 


+61 10 58.998(12) 


3.1(2) 


2.9(3) 


150(110) 


<250 


150(24) 






A- 


-array 










W 33A Ql 


18 11 44.2104(4) 


-17 52 57.804(9) 


0.80(17) 


1.7(5) 


87(57) 


<67 


41(20) 


W 33A Q2 


18 11 44.2510(3) 


-17 52 58.241(6) 


0.79(18) 


0.62(27) 


<49 


<12 


a 


W 33A Q3 


18 11 44.1769(3) 


-17 52 58.972(5) 


0.76(18) 


0.53(25) 


<37 


a 


a 


AFGL 2591 


20 27 35.9520(3) 


+40 01 14.682(4) 


0.95(16) 


1.9(5) 


50(15) 


29(23) 


33(39) 


NGC 7538 IRS9 


23 11 52.8671(3) 


+61 10 59.017(4) 


1.1(2) 


1.1(3) 


40(17) 


a 


173(17) 



a : No solution was found for this parameter 



sources are only resolved in one direction. In the few cases 
where no meaningful estimate of the minor axis could be 
obtained, we use the value for the major axis as upper 
limit. Using the distances from ^ the A-array sizes cor- 
respond to linear radii of 150 AU for W 33A, 20 AU for 
AFGL 2591, and 55 AU for NGC 7538 IRS9. 

The sizes found from the A-array data are not con- 
sistent with the C-array estimates, but factors of 5-7 
smaller. This discrepancy suggests that the sources do not 
have a uniform brightness, but rather a centrally peaked 
brightness distribution. Each array picks up the 'tip of the 
iceberg' at its own level of brightness sensitivity. 

The derived axis ratios (or limits) of <0.5 suggest 
a nonspherical, flattened geometry. The position angles 
found from the A-array images are consistent with those 
found from the C-array images within their 2a errors. 
Combining the results from both array configurations, the 
orientation on the sky is PA«20° for W 33A Ql, PA«73° 
for AFGL 2591 and PAwl65° for NGC 7538 IRS9. 

3.2. Flux densities and brightness temperatures 

Columns 4 and 5 of Table Q give the peak intensities and 
the integrated flux densities of the sources. The flux den- 
sities determined from the A-array data are «two times 
lower than those derived from the C-array data. Since the 
phase calibrators are not systematically brighter in either 
configuration, we rule out a systematic calibration offset 
between the two data sets. 

Sources can be weaker in more extended configurations 
because they are resolved, or because of increased phase 
noise on longer baselines. This 'atmospheric decorrelation' 
decreases the visibility amplitudes by a factor of e -0 */ 2 , 
where o$ is the phase noise (in radians). From the phase 
noise in the calibrator data, we estimate a decorrelation 
by «5% in C-array and by «10% in A-array. Since the 
discrepancy is much larger than these factors, we conclude 
that decorrelation is unimportant here. Combined with 
the difference in observed sizes (S I3.1|) . this discrepancy in 



flux density is additional evidence for a centrally peaked 
brightness distribution for these sources. 

A first clue to the nature of the emission is the conver- 
sion of the measured flux densities and deconvolved sizes 
to brightness temperatures, using the Rayleigh- Jeans law. 
The total flux densities measured with the C-array cor- 
respond to Tb— 20 - 50 K, while the values measured 
in A-array indicate Tb=190 - 850 K. Most of these num- 
bers are underestimates of the true brightness because the 
emission is not fully resolved, so we adopt Tg — 500 - 
1000 K as representative range for these sources. 

If due to ionized gas, which has an intrinsic tempera- 
ture of 10,000 K, the emission must be beam diluted, or 
optically thin, or both. However, the cm-wave spectra of 
W 33 A and AFGL 2591 (§ E3J are too steep to be due 
to optically thin free- free emission. Therefore, if the ob- 
served emission is from ionized gas, the filling factor must 
be ~0.05. Using the source sizes derived in 5 13.11 this fill- 
ing factor implies a linear radius of «30 AU for W 33A, 
«5 AU for AFGL 2591 and «10 AU for NGC 7538 IRS 9. 

If the 43 GHz emission is due to dust, the maximum 
physical temperature is «1500 K, where dust sublimates. 
Therefore the emission must be close to optically thick 
and close to resolved. The mm-wave spectra (S I3.3fl indi- 
cate optically thick emission in W 33 A and AFGL 2591, 
but optically thin emission in NGC 7538 IRS9. The fill- 
ing factor must be at least 1/3 and is probably close to 
unity, especially if the emitting dust is cooler than 1500 K. 
The corresponding physical sizes assuming single conden- 
sations are 90 - 150 AU for W 33A, 5 - 20 AU for AFGL 
2591, and 10 - 30 AU for NGC 7538 IRS 9. 

3.3. Radio spectra 

Table contains radio flux densities for our three sources 
from the literature. Only data between frequencies of 3- 
300 GHz and in beams < 3" have been included. Figure^] 
plots the data from Tables and [3 We use the total flux 
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densities measured in C-array, which are most comparable 
to the literature data. 

The broad-band spectrum of W 33A at frequencies up 
to 43 GHz can be characterized by a spectral index 7, 
defined as 

S v oc v 1 

with 7=1.2±0.4, while the spectrum steepens to 
7=2.3±0.2 at higher frequencies. Similarly, for AFGL 
2591, 7=0.9±0.4 at v < 43 GHz, while 7=2.5±0.8 at v > 
43 GHz. The spectrum of NGC 7538 IRS9 is even steeper, 
7=3.0±0.5, at v > 43 GHz, and inverted (7=-0.5±0.4) at 
centimeter wavelengths. The uncertainties in these spec- 
tral indices are dominated by the differences in values 
derived at different frequencies, rather than by the un- 
certainties of the individual flux densities (estimated as 
10% for the VLA, 30% for OVRO at 86 GHz and 50% at 
230 GHz). The spectral indices may be affected by vari- 
ability, since the mea surements were not simult aneous. In 
the case of W3 IRS5. IVan der Tak et all ll2005h found ra- 
dio variability by factors up to 3. 

The spectral indices of W 33A and AFGL 2591 in- 
dicate that the emission at z^<43 GHz is dominated by 
free-free emission, while the higher-frequency points are 
mostly due to dust emission. Similar conclusions hold for 
AFGL 2 136 and Orion-L which indee d have similar radio 
spectra l|Menten fc van der Tak1l2004|) . 

The spectrum of NGC 7538 IRS9 differs from that of 
the other sources in two respects. First, the mm- wave spec- 
trum is significantly steeper than 7=2, indicating that the 
dust emission in NGC 7538 IRS9 is optically thin. Second, 
NGC 7538 IRS9 has a 'flatfish' spectrum between 4.9 and 
8.5 GHz, whose spectral index is consistent with 0, but 
poorly determined, due to the limited sign al-to-noise ra- 
tio in particular of the 8.5 GHz emission l|Sandell et alJ 
l2005h . The latter authors interpret this as free-free emis- 
sion, which, given the spectral behaviour would have to be 
optically thin, which is unlikely for such a compact source. 
The sizes they determine for the radio source, 2" x 1" 
at 8.46 GHz and 0.7" x 0.1" at 4.86 GHz, also are very 
uncertain, and should probably at best be considered as 
upper limits (consistent with our 43 GHz data). We note 
however, that within the uncertainties, NGC 7538 IRS9's 
radio spectrum could be due to non-thermal emission, con- 
sistent with synchrotron emission from a jet, as seen be- 
fore in the 'Turner- Welch' object W3(H 2 0) bv lReid et all 
l)l995(l . An origin of the radio emission in a jet is consis- 
tent with the source's morphology, which is unresolved in 
one dimension. Its orientation, almost north -south, is sim- 
ilar t o that of the high- velocity CO outflow ijSandell et al.l 
2005j), although we note that outflows around NGC 7538 
IRS9 observed in other tracers (H2 and HCO + ) present a 
complex picture. 

The centimeter-wave spectral indices of W 33A and 
AFGL 2591 are intermediate between the values expected 
for optically thin (7=-0. 1) and optically thick (7=2) H II 
regions with a uniform temperature and density. An in- 
termediate optical depth would lead to a 'bent' spectrum 



Table 2. Overview of literature data. 



Source 


V 

(GHz) 


s v 

(mjy) 


Reference 


W 33A 


8.4 


0.8 


4 




15 


1.9 


4 




86 


24 


5 




233 


190 


5 


AFGL 2591 


4.9 


0.4 


1 




8.4 


0.5 


2 




22 


1.6 


6 




87 


30 


3 




226 


151 


3 


NGC 7538 IRS9 


4.9 


1.0 


7 




8.5 


0.8 


7 




15 


<0.5 


4 




107 


43 


5 



(1) ICampbelll Il984l (21 iTofani et alJ I l995 | 
' ' 119991 (4) iRengaraian fc Hcl 19961 
........ . I2OOOI (6) iTrinidad et al] 12003 (7) 

ISandell et al.ll2005l 



References: 

(3) IVan der Tak et alJ 
(5) IVan der Tak et alJ 



which is not observed. However, the observed values of 
7 ss 1 are expected in the case of H II regions with den- 
sity gradients. The most plausible gradients are the cases 
of an ionized wind and an infalling envelope, which lead to 



and s 



-1.5 



density structure, respectively. The 



expected spectral indices of 7=0.6 and 7=0.1 are some- 
what below the observed values. 

In the case of a power-law distribution of the electron 
density with radius, n e oc r~ 9 , a spectra l index of 7 = 
(2g-3.1)/(g-0.5) is expected l|01nonll975|) . Therefore the 
observed value 7 ~ 1 would indicate q — 2.6. However, it 
is not clear which process would create and sustain such a 
steep density distribution. The sound speed of «10 km s _1 
implies a crossing time of «50 yr for our sources, so the 
process would have to act continuously. More plausible 
models will be discussed in § 14.31 



4. Discussion 

4.1. Stellar ionization rates 

The flux densities S v in Table^can be combined with the 
distances d from §Hto estimate the flux of Lyman contin- 
uum photons Nl emitted by the central star fe.g.. lKrugell 
1200.1) . In the optically thick limit, the 'Stromgren' radius 
rs of the H II region follows from blackbody emission at 
an assumed electron temperature of T e =10 4 K: 

7T 

S v =B v (T e )-£ — 

The derived radii of 33.5 AU for W 33A, 7.3 AU for 
AFGL 2591, and 19.3 AU for NGC 7538 IRS9 are consis- 
tent with the observational limits (§ 13.1(1 . To have a free- 
free optical depth of unity, the emission measure must be 
at least 8.4 xlO 9 cm~ 6 pc. These radii and emission mea- 
sures together imply electron densities of n e =5xl0 6 to 
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Fig. 4. Filled squares: Small-scale (<3") radio (3- 
300 GHz) spectra of W 33A (top), AFGL 2591 (mid- 
dle) and NGC 7538 IRS9 (bottom). The 'ionized accretion 
flow' model from § 14.31 (dotted lines) is seen to reproduce 
the emission up to ^50 GHz to within a factor of two. Flux 
densities at higher frequencies >50 GHz are dominated by 
dust emission, which the model does not consider. 



lxlO 7 cm 2 . The ionizing photon flux then equals the 
recombination rate: 



4tt o 2 

N L = -J r S n e aB 



where cub is the 'case B' recombin ation coefficient of 
1 l)Osterbrocklll989|) . The values of N L 
1 for AFGL 2591, 
Values estimated 



are3.5xl0 45 s" 1 for W33A, 1.7xl0 44 s" 
and 1.2xl0 45 s" 1 for NGC 7538 IRS9. 
using the optically thin limit, 



Nr. 



7.54 x 10 



46 ' 



»_»/___ Y (JL.Y ' 45 ' " 

Jy Vkpc / V 104k / 
are within a factor of 2 from the thick estimates. 



(gHz) 



4.2. Dust accretion rates 

Ionizing photon fluxes can also be estimated from the in- 
frared luminosities, assuming that one single star dom- 
inates the Lyman continu um flux. We use the Gen eva 
stellar structure models l)Maeder fc MevneH Il989h to 
convert luminosities into stellar masses and effective 
temperatures, and t he st ellar atmosphere models by 
ISchaerer fc de Koterl Jl997h to estimate Lyman contin- 
uum fluxes. The results are M=28 M Q , T off =40,000 K and 



A L =8.5xl0 48 s- 1 for W 33A, M=16 M Q , T cff =33,000 K 

^48 



and Af L =1.0xl0 4 " s" 
T eff =35,000Kand A L = 



for AFGL 2591, and M=20 M , 
: 2.2xl0 48 s- 1 for NGC 7538 IRS9. 
These ionizing photon fluxes are factors of 2000-6000 
higher than the estimates from the radio continuum emis- 
sion (S I4.1|I . Dust inside the H II region only accounts for 
factors of 2 - 3 absorption. One possibili ty is confinemen t 
by the ram pressure of accreting dust (cf . I Walmslevll 1 9 9 5h . 
The required accretion rate follows from equating the ram 
pressure to the thermal pressure of the H II region: 

Mvft = 4:TTr 2 s n e k B T e 

where the free-fall velocity is given by 



vg = ^GM^/rs 

which equals 42 km s" 1 for M*=20 M Q . With r s =20 AU 
and n e =10 7 cm -3 , the estimated accretion rate is 6xl0 -8 
M yr~ 4 . This value is a lower limit because radiation 
pressure may slow down the dust from the free-fall speed, 
and because a stellar wind may add to the internal pres- 



4.3. Ionized accretion flows 

Recent work bv lKetol ll2002b|) shows that stellar gravity 
prevents the hydrodynamical expansion of H II regions 
inside a 'gravitational radius' 

r g = GM/2c 2 s 

where c s is the isothermal sound speed of ionized gas 
of wlO km s _1 . This radius lies at 124 AU for W 33A, 
71 AU for AFGL 2591, and 89 AU for NGC 7538 IRS9. 
These values are comparable to the observational limits 

(§EJ. 

In Keto's model, gas accreting onto the star is in molec- 
ular form at large radii, and becomes ionized at r = r g . 
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The density structure of the molecular and ionized regions 
is given by the n oc r -15 free-fall profile. The molecular 
envelopes of our s ources have similar or s lightly flatter 
density structures l|Van der Tak et aDl2000|) . 

The calculation of the radio spectru m of g ravita - 
tionally bound H II regions is given by IKetol l|2003j) . 
Appendix 1X1 summarizes the derivation, with a few errors 
in the formulae corrected. The dotted lines in Figure^lare 
flux densities calculated using these corrected formulae. 

The only free parameter in the calculation is no, the 
density at the molecular sonic point r m . Adopting T=30 K 
as representative for the molecular gas, the location of 
its sonic point only depends on the mass of the central 
star. We have fitted the parameter n to the observed 
(unweighted) flux densities of our sources at frequencies 
up to 43 GHz. 

For W 33A, we find n =3xl0 5 cm" 3 at r m =0.49 pc, 
for AFGL 2591, we find n =8xl0 4 cm- 3 at r m =0.28 pc, 
and for NGC 7538 IRS9, we find n =2xl0 5 cm" 3 at 
r m =0.35 pc. These densities agree to factors of a few with 
the values derived fr om maps of sub-millimete r continuum 
and molecular lines l|Van der Tak et ai~l |2000). 

Figure0]shows that the model reproduces the observed 
flux densities to a factor of ~2, which is the expected 
uncertainty due to variability (in particular, variations in 
the accretion rate). We conclude that ionized accretion 
flows describe our data well, but that stronger tests are 
necessary to prove the model beyond doubt. 

4.4. Dust disks 

As an alternati v e mod el, we consider the suggestion by 
IPreibisch et al.l l)2003|) that the near-infrared emission 
from AFGL 2591 is the inner rim of a thick disk or enve- 
lope at the dust sublimation radius. These authors resolve 
the near-infrared emission and derive a diameter of 40 AU, 
which is in good agreement with our measured 43 GHz 
size, assuming a beam filling factor of unity. Since our 
measured 43 GHz brightness is 850 K rather than 1500 K, 
this interpretation requires either the optical depth to be 
w0.5 at 43 GHz, or the dust temperature to be «850 K, 
about half the sublimation temperature. 

The luminosity of a T=1500 K, 40 AU diameter black 
body sphere is 80,000 L , which exceeds the observed 
value of 20,000 L Q (§^l. Therefore, if the object has this 
temperature and size, it can fill at most 1/4 of the sky, and 
is likely a disk. However, if T=850 K, the spherical black 
body luminosity is only 9000 Lq , within the observational 
limits. 

The requirement of unit optical depth at 43 GHz con- 
strains the dust colu mn density. Extrapolat i ng Mo del 5 of 
the dust models by lOssenkopf fc Hennind l)l994|) , which 
provides a g ood fit to the mid- to far- infrared spectrum of 
AFGL 2591 l|Van der Tak et al.ll999|) . the 43 GHz opacity 
is 0.06 cm 2 g _1 . The 2.2 /mi opacity is a factor of 80,000 
higher in this dust model, so that t(43 GHz)=1 implies 
Ay=890,000 or AT H =1.8x 10 27 cm" 2 . This column density 



is orders of magnitude higher than m easured for the en- 
velopes of embedded high-mass sta rs l|Van der Tak et all 
l2000UHatchell fc van der Takll2003h . but may be reached 
by massive circumstellar disks. 

Another check of the disk model lies in the 2.2 /im 
flux density of 2.8 Jy observed toward AFGL 2591 
(e.g.. lAitken et ahlll988t) . The compact component in the 
speckle observations contributes «60% of the 2.2 /Ltm flux 
density of AFGL 2591 in a 12" field (T. Preibisch, priv. 
comm.). If the emission is a black body at T=1500 K 
with a radius r=20 AU, matching the observed bright- 
ness requires a foreground 2.2 /jm extinction of 25 mag- 
nitudes, or 20 mag if T=850 K. The corresponding visual 
extinctions of 200-250 imply N H =2.0-2. 5 x 10 23 cm" 2 , 
which agrees very well with the N(H2 ) = 2.3xl0 23 cm" 2 
meas ured in sub-millimeter emission l|Van der Tak et alJ 
1999 D. We conc l ude th at the dust disk seen at 2.2 fim by 
IPreibisch et alJ (|2Q03) may also be responsible for our ob- 
served 43 GHz emission. The evidence for this interpreta- 
tion would be strengthened by measuring the orientation 
of the near-infrared emission, which the current speckle 
data do not constrain (T. Preibisch, priv. comm.). 

4.5. Relation with large-scale outflows 

The deconvolutions in § 13.11 indicate flattened shapes for 
the radio sources, with axis ratios of «3. One central pre- 
diction of disk accretion models of low-mass star forma- 
tion is that the outflow axis is perpendicular to the disk 
plane. To see if such models may also work for high-mass 
star formation, this section compares the orientation of 
the 43 GHz emission with that of the large-scale outflows 
observed in our sources. 

It is difficult to constrain the outflow orientation of 
such objects from the commonly observed CO 1-0 line. 
The outflows appear poorly collimated in single-dish ob- 
servations, and interferometric observations tend to re- 
solve out most of the emission. For example, in AFGL 
2591, one of the earliest known outflow sources, high ve- 
locity CO in the 1-0 line is spread ove r a large area (~5 7 
or ^1 pc) and is poorly collimated l|Lada et alJ Il984j) . 
Better constraints come fr om higher— J CO l ines or other 
molecules such as HCO + (|Benther et alJl2002^ . 

In the c ase of AFGL 2591, CO 3-2 and HCO+ 4-3 
mapping by iHasegawa fc Mitchelll (H995) shows an out- 
flow of size 90x20", embedded in an arcminute-scale flow 
seen before in lower- J CO lines. The East-West orien- 
tation agrees with the positions of Herbig-Haro objects 
and spots of shock-excited H? emission. For NGC 7538 
IRS9, maps of the CO 2-1 line l|Mitchell et alJll99lh show 
an outflow with a North-South orientation (PA~160°). 
Very high velocities (up to 110 km s _1 ) are seen in 
the 3-2 line, but this emi ssion is spatially unresolved 
l|Mitchell k. Hasegawalll99l|) . For W 33A, outflow activ- 
ity is demonstrated by the positions and velocities of OH 
and H2O masers (e.g., iFish et al1l2003l) and by 'wings' 
on the sub-millimeter emission line profiles of CS, HCO + 
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and other molecules (van der Tak et al. 2000, 2003). 
Unfortunately, there is no CO or HCO + map to give us 
the outflow orientation for this source. In fact, the out- 
flow may not be bipolar, as OH masers usually appear in 
quasi-spherical outflows. Observations of H 2 masers with 
0.5 mas resolution have revealed spherical outfl ows from 
sever al other young high-mass stars fe.g.. iTorrelles et all 
12003ft . 

The 43 GHz emission in both AFGL 2591 and NGC 
7538 IRS9 has a similar orientation on the sky as their 
large-scale outflows. This alignment argues against an ori- 
gin of the 43 GHz emission in a disk, and suggests an 
origin in a jet. However, thermal emission from an ionized 
jet su ch as seen in lower- mass objects fe.g.. lTrinidad et all 
l2004h is inconsistent with the measured spectral indices 
(§ 13.31) . Perhaps cool material entrained by the jet con- 
tains sufficient dust to produce the observed spectrum. 



4.6. Photo-evaporating disks 

Another potential source of our observed 43 GHz emission 
is free-free emission from an ionized wind from a circum- 
stellar disk that is being illuminated by the Lyman con- 
tinuum of a hot central star. Such photo- evaporating disks 
have b een described semi-analytically by Hollcnb ach et alJ 
(1994) while numerical calculations were performed by 
iRichhng fc Yorkellll997ft . In these models, evaporation oc- 
curs in the outer disk, where the radiative force of the 
star is stronger than its gravitational force. In our case of 
a ^20 M Q star, this radius occurs at ~100 AU, which is 
comparable to our measured source size. 

The models distinguish 'strong' and 'weak' stellar 
winds, depending on whether the ram pressure of the wind 
is comparable to the thermal pressure of the disk atmo- 
sphere or exceeds it. Wind speeds up to several 100 km s _1 
have been measured in mid-infrared CO absorption and 
H emission lines toward AFGL 2591 an d NGC 7538 IRS9 
I Mitchell et "aTl ll 99ll lBu"rm" et al.l ll995l) and may also ap- 
ply to W 33A. For a wind speed of 500 km s _1 , a stellar 



mass of 20 M Q and a Lyman continuum flux of 3 x 10 s 
f§ I4.2fl . the critical mass loss rate is 2xl0 -7 M Q yr _1 . 
The mass loss rates estimated from CO 2-1 are ~10 -4 
M^yr -1 , so the sources are clearly in the strong wind 
limit. 

The predicted spectral shape of circumstellar disks 
that arc being photo-evaporated by strong stellar winds 
is that of an optically thin H II region (7=— 0.1) at short 
wavelengths, a thick H II region (7=2) at long wave- 
lengths, and a 'wind' spectrum (7=0.6) at intermediate 
wavelengths. For the high mass loss rates of our sources, 
the transition to optically thick emission lies at unob- 
servably long wavelengths. The transition from 7=0.6 to 
7=— 0.1 is expected around 0.5 cm wavelength. However, 
the observed spectra have quite different shapes, so that 
the photo-evaporating disk model does not seem applica- 
ble to our sources. 



4. 7. Free electron absorption 

As a last option, we consider free-free absorption by H as 
source of the observed 43 GHz emission. This mechanism 
was originally prop osed for the radio em ission of Mira- 
type variable stars llReid & Mentenlll997 ) and may play 
a role in Orion-KL IjBeuther et alJ 12004 K The required 
free electrons come from photo-ionization of Na and K 
atoms, so that this mechanism only operates inside the 
dust sublimation radius. 

Our observed brightness temperatures of ~1500 K 
are consistent with this model, and combining the radii 
of 
TO 



10 AU with A^(H2)~10 27 cm" 2 gives a density of 



iu 13 c m" 3 , high enough to p rovide unit free electron 
opacity ijReid k, MenterJ Il997[) . The broad-band radio 
spectra of W 33 A and AFGL 2591 are not inconsistent 
with the predicted 7=1.86 behaviour if the flux densities 
are uncertain by a factor of ~2. For NGC 7538 IRS9, free 
electron absorption may explain the 15 - 107 GHz spec- 
trum, while the lower-frequency data must have a different 
origin. In any case, stronger tests of this model will come 
from measurements of the density on ~T0 AU scales. 



5. Conclusions 

We have detected weak, compact radio emission from 
three high-mass protostars which are deeply embedded in 
molecular envelopes. In two cases, the emission is prob- 
ably due to bremsstrahlung from gravitationally bound 
H II regions. The derived densities are in good agreement 
with values found for the large-scale molecular envelope. 
Similar results were foun d for another object, W3 IRS 5 
ilVan der Tak et al.ll2005|L and additional support for ac- 
cretion onto high-mass pro tostars comes from H?Q mase r 
mapping of AFGL 2136 (|Menten & van der Takl l2004h . 
These findings support a picture of high-mass star for- 
mation where the central star builds up mass by accre- 
tion from a molecular envelope. At small radii, the stel- 
lar Lyman continuum ionizes the accretion flow, but does 
not stop it. As the star gains mass, its surface tempera- 
ture and ionizing flux increase. The circumstellar ionized 
region grows quasistatically until its radius reaches the 
critical value of GM/2c? s . Stellar gravity is then unable to 
confine it, so that the H II region rapidly expands and 
halts the accretion. 

However, this interpretation is not unique and may not 
be complete. For AFGL 2591, the data are also consis- 
tent with dust emission from a disk seen previously in the 
near-infrared. Furthermore, NGC 7538 IRS9 appears to 
have a non-thermal cm-wave spectrum, such as produced 
by synchrotron emission from a jet. The mm-wave spec- 
trum of this source is steeper than that of the other two, 
suggesting optically thin dust where the others are thick. 
In the case of W 33A, three radio continuum sources are 
seen at 3000 - 5000 AU separations. This 'clustering' of 
radio sources resembles the W3 IRS5 region (<i=2.0 kpc, 
L=1.4xl0 5 L©), where seven sources are found a t abou t 
half the above separations l|Van der Tak et al ] l2005h . 
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Follow-up VLA observations would be useful to find out if 
some of the 43 GHz sources in W 33A are transient, and 
mid-infrared imaging would reveal which of the three are 
self-luminous . Together with our previous 'case studies' of 
AFGL 2136 llMenten fe van der Takll2004|) and W3 IRS5 
llVan der Tak et alJl2005lT the current observations sug- 
gest that the formation of stars and stellar clusters with 
luminosities up to ~10 5 L Q proceeds through accretion 
and produces collimated outflows as in the case of solar- 
type stars. Such region s may well contain small star clus- 
ters llTesti et al . 1999) and clusters of non-thermal radio 
sources l|G6mez et al"[l2002j) . The main difference n the 
formation of ~10 4 L Q stars with that of solar-type stars 
is that their accretion flows become ionized at small radii. 
In contrast, regions of ~10 5 L Q appear to produce clusters 
of H II regions and essentially uncollimated outflows, pos- 
sibly as the result of mergers of protostars or pre-stellar 
cores. 

In the future, the ionized accretion flow model should 
be tested by measuring the density profiles of hypercom- 
pact H II regions directly through radio continuum ob- 
servations at ~10 AU resolution (10 mas at 1 kpc) and 
~0.1 mjy sensitivity. The duration of the ionized accre- 
tion flow phase should be measured by a radio continuum 
survey of high-mass protostars at a resolution and sen- 
sitivity comparable to this work. Observations of radio 
recombinatio n lines at high resolution exist for one source 
(|Ketoll2002a|) and should be carried out for more. These 
efforts probably exceed the capabilities of the VLA, but 
will be feasible with the e-VLA (after 2006) and ALMA 
(after 2008). 

Acknowledgements. The authors thank Friedrich Wyrowski, 
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If the emission is unresolved, the total flux density is 
S v = 2fc B (T t hick + T thin )Oj/ 2 /c 2 

where 

fl = n9l 

is the solid angle of the H II region. 



Appendix A: Flux density of an ionized accretion 
flow 

An H II region with a density gradient consists of an op- 
tically thick core and an optically thin envelope. The an- 
gular radius of the optically thick core 9\ is related to the 
angular radius of the H II region 8 by 



with 



and 



B = A FF n 2 T e 9 2 — 
pc 



Y FF = 8.235 x 10- 2 T- 35 (^-) 



-2.1 



,GHz; 

The radio emission from a gravitationally bound H II 
region is the sum of the emissions from the optically thick 
core and the optically thin envelope, weighted by their 
respective areas on the sky. For the optically thick core, 
the brightness temperature is 

q2 

^thick — T e —77 

while for the optically thin envelope, 

d ( 6 
Tthin = -AA FF nl— u + 9 /2 In 

pc V #0 

with 



